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We consider the evolution of a limited Quark-Gluon Plasma (QGP)drop in the 
framework of relativistic hydrodynamics. In the presence of the boundary, the ex- 
panding and cooling of QGP drop may appear as oscillating processes, which realize 
via multiple traveling of rarefaction and compression waves through the plasma and 
leads to increasing of plasma life-time. Within this picture, deflagration seems to be 
the main possible way of hadronization. The slow cooling mechanism also leads to 
a significant enhancement of photon yield in comparison with Bjorken scenario. 

1. Introduction 

When discussing the signals of a quark-gluon plasma produced in heavy 
ion collisions, particularly the ones connected with thermal photons or 
lepton pairs [jlj], it is necessary to have an adequate underlying picture 
of the space-time evolution of plasma. Moreover, the very characteristic 
details of the evolution of plasma may be considered as signals by their 
own 

In principle, it is generally accepted that relativistic hydrodynamics is 
able to provide the wanted description, but the exact solution of the math- 
ematical problem is rather difficult. The approach, which is most widely 
exploited in this context, if not to say the only one, is known as Bjorken 
model This model is based on the assumption of a space-time scaling, 
which means that all processes in a hot strong-interacting system do only 
depend on the local proper time r. The assumption of space-time scaling 
leads to major simplifications and allows to obtain an analytical solution 
of the hydrodynamical equations. At a further step, this scaling solution 
is often used as a ground for more sophisticated models ^ where the 
transverse motion of matter is taken into account as a correction to its 
scaling-like longitudinal expansion. 

However, this model is only an idealization, which ignores some impor- 
tant features of the real world. For example, it is the fact that the QGP 
drop is, anyway, limited, and so there must exist a boundary, which sepa- 
rates two phases with different equations of state. 

Our own study is based upon a direct numerical analysis of hydrodynam- 
ical equations together with the boundary conditions between plasma and 



vacuum. In the present work we have considered the space-time evolution 
of a spherically-symmetric QGP drop and have tested various sets of initial 
conditions. 

A similar approach, also applied to a spherically symmetric QGP drop, 
has been proposed in ref. []7|. The author of []7| mainly focuses his atten- 
tion on the effects of pion-pion Bose-Einstein correlations and compares 
his predictions with Bjorken Model. This is well complementary to our 
present consideration, where the emission of thermal photons is revised. 
In addition, we include in our analysis various ways of hadronization, both 
equilibrium and nonequilibrium, whereas the paper [)7| deals with only one 
type of phase transitions, the deflagration. 

Having the a priori space-time scaling given up, we get a new, quali- 
tatively different view of the evolution of plasma. Among the noticeable 
consequences of our analysis is the fact that supersonic phase transitions 
from QGP to hadrons are impossible, and this may strongly affect the even- 
tual proportion between different particles, which appear in the course of 
hadronization [p, |^ . Another important observation is the enhanced yield 
of thermal photons, which exceeds the traditional estimations by about 
one order of magnitude. Our results remain valid in a wide range of model 
parameters. 

The outline of the paper is the following. In sect. 2 we present the basic 
equations of hydrodynamics and the equation of state of the deconfined 
matter. In sect. 3 we briefly discuss the numerical scheme to integrate the 
dynamical equations with boundary conditions. In sect. 4 we analyze the 
evolution of a QGP drop as it appears in our approach and point out some 
essential features of the dynamics. In sect. 5 we calculate the production 
of thermal photons. Our main findings and conclusions are summarized in 
sect. 6. 



2. Equations of motion and equations of state 

To describe the evolution of matter created in a relativistic collision of 
heavy nuclei one commonly uses the hydrodynamical approach [|TD|, [TT|. 
The equations of motion represent the momentum and energy conserva- 
tion. In the absence of dissipative processes (diffusion, viscosity, heat con- 
ductivity) they read 

5^r^- = 0, (1) 
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where the energy-momentum tensor 



T^^ = {£+p) u^'u" - pgt'". (2) 

Here e is the energy density, p is the pressure, is the 4-velocity, and g^'^ 
is the metric tensor. 

For the spherically symmetric case considered here it is convenient to 
rewrite the above equations (1),(2) in an explicitly symmetric form: 

,1 n / 2x./ NX n -('^^^-P] ff . _ ({w-f'^ - p)Vr+pVr 

dtu + ^drir i{u)) = 0,u^ \ , f u) = 



where w = £ + p is the enthalpy density, is the radial component of the 
velocity, and 7 = (1 — v'^)~^^'^. 

The dynamics of a QGP drop is now completely defined by the initial 
conditions, the boundary conditions and the equation of state, which re- 
lates p and £. 

We adopt the MIT-motivated equation of state of plasma []T^ and get 
for massless quarks and gluons 

p = aT^- B, 
£ = SaT"^ + 5, 

where T is the temperature and 5 is a phenomenological constant corre- 
sponding to the energy density and the pressure of the QCD vacuum. The 
parameter a defines the number of degrees of freedom in QGP and is taken 
a = 7rV90(2 x8 + 2x2x3xiVjx|)« 237rV45, where Nf = 2.5 for 
taking into account u- , d-, s-quarks. 

Since the QGP drop is assumed to be of limited size it is surrounded by 
the hadronic vacuum with the equation of state 

^vac Pvac 0- (5) 

We will recall this equation when formulating the boundary conditions. 



3. Numerically integrating the hydrodynamical equations 

Numerical algorithms aimed at the solution of the hydrodynamical prob- 
lems are being developed for many years. Among the various approaches 
described in the literature we prefer the nonlinear monotonic transport 
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scheme FCT (Flux-Corrected Transport) [jT^, [TJ, [T5j] developed by J.P.Boris 
and D.L.Book in the framework of Eulerian grid methods. It has been 
shown to apply to a wide class of physical systems and is found quite 
accurate and reliable. 
Let us remind the guideline of this method. 

Consider one-dimensional continuity equation for the mass density flux: 

dtp + d,{pv) = 0. (6) 

Let denote the density in the i^^ point of the grid at a given time. An 
explicit scheme to calculate the density pi at the next time step may be 
presented as follows: 

Pi = Z^? - Q+i/2(^°+i + + Q-i/2(/?°-i + + 

^l+l/2{P^^+l - P^) + ^l-l/2{P^^-l " P^) ■ (7) 

This form observes the mass conservation. Here the quantities with half- 
integer indices ? ± 1/2 correspond to arithmetic-mean values of the ones 
taken in the points i and ? ± 1, and €^±1/2 — 0-/'^) Vi±ii2{5t/ 5x) with 5t 
being the time step and 5x the space constant of the grid. The terms 
with V stand for an artificial diffusion and are introduced to guarantee the 
positivity of mass density, which demands 

^i±i > kiiil for all i (8) 

To compensate this overstrong diffusion one applies an antidiffusion cor- 
rection: 

pI = Pi- Pi+i/2{pi+i - Pi) - Pi-i/2{pi-i - Pi), (9) 
with Pi±i/2 being positive coefficients. This step removes extra smoothing 
but no longer guarantees the positivity of the newly computed densities 
pj. To restore the positivity and the stability of the eventual solution the 
antidiffusion fluxes fi±i/2 = Pi±i/2{Pi±i — Pi) need to be corrected. In the 
simplest form one can take Pi±i/2 = ^i±i/2—\^i±i/2\ - Eliminating the residual 
diffusion [v — p) is only possible with nonlinear corrections adjusted to the 
local behaviour of the solution. 

The fundamental requirement of the FCT method is that the antidiffusion 
step must neither generate new extrema in the solution nor enhance the 
ones already existing. This is observed with the following prescription: 

1 ~ fcor fcor 

Pi — Pi Ji+l/2 Ji-1/2^ 

fi±i/2 = 'S'maxjO, mm[S{p,±2- p^±l), \f^±l/2\, S {p, - p,^i)]} , (10) 
l^l = 1, sign^ = sign(p,±i - p,). 



4 



The condition ( pX)|) is constructed in such a way that it takes into account 
all possible combinations of the signs of local gradients of the physical 
solution. A detailed description of various transport, antidiffusion and 
correction schemes may be found in [|T5|. Like it is done in the original 
papers [0, [T5|, we also use the time-splitting procedure to calculate the 
velocities. 

An essential point in our analysis is the presence of the moving phase 
boundary, which greatly complicates the consideration. In general, La- 
grangian grids are better suited for this kind of problems. Since the grid 
cells are made moving together with the matter, one can identify the in- 
tercell boundary with the phase boundary, thus avoiding the ambiguous 
situation of mixing two different equations of state in one cell. To fit in 
with the FCT method we have adopted a combined scheme. 

All the internal cells of our grid are of Eulerian type and have definite 
unchangeable size (it will be refered to as "standard size"). On the con- 
trary, the cell adjacent to the phase boundary is of Lagrangian type, with 
its outer face moving together with the boundary. Under the evolution of 
the QGP drop the size of this cell changes continuously. If, due to the 
expansion of the drop, the length of the surface cell becomes larger than 
1.5 standard lengths, the cell splits into two new cells. The one, which 
does not touch the phase boundary, gets the full standard size. It now 
belongs to the Eulerian part of the grid. The other cell, which is adjacent 
to the phase boundary, takes the rest of the total length and continues its 
life as Lagrangian cell. In the opposite case, if the length of the surface 
cell reduces to a quantity less than 0.5 standard lengths, it fuses with the 
neighboring Eulerian cell, thus forming a Lagrangian cell of extended size. 

All the internal Eulerian cells are treated according to the standard FCT 
algorithm. The parameters of the boundary Lagrangian cell are found from 
the momentum and the energy conservation. The incoming and outcoming 
fluxes of the momentum and the energy for the Lagrangian cell are known 
from the parameters of the neighboring Eulerian cells on the inner side 
and from the equation of state of vacuum (g) on the outer side. Given 
these fluxes, we calculate the new values of the momentum and the energy 
densities. Then we extract from these quantities the velocity of plasma. 
This velocity is assumed to be the velocity of the boundary of the QGP 
drop. Given the velocity, we calculate the new size of the Lagrangian cell 
and recalculate the momentum and the energy densities according to the 
changed size. 
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The calculated values of the densities and the volume of the Lagrangian 
cell may then be corrected in such a way that to achieve an exact conser- 
vation of the energy. This means 

N N~l N-1 

£nVn = E ^i^i - E ^i^t = E^t=o) - E ^i^i^ (11) 

i—l i—l i=l 

where £'(f=o) is the initial total energy of the system. This equality may be 
used either to correct the volume V/v of the Lagrangian cell (if to think the 
energy density e^r is known correctly) or to correct the energy density (if 
to think the volume is known) or to make a weighted correction for both 
these quantities. The latter case is realized in our paper. This scheme was 
tested in 1-dimensional case with very simple Lagrangian method [ITO . 



4. The dynamics of a QGP drop 

In the present section we consider the dynamical properties of the evo- 
lution of a QGP drop. The initial conditions for this problem must be 
derived from the dynamics of nucleus- nucleus collisions. In fact it leaves 
a large piece of freedom. We do not pretend to give precise and detailed 
predictions on the behaviour of plasma in real experimental conditions. 
Our purpose is rather to give a qualitative description, which reveals the 
characteristic features of a new mechanism of expansion and cooling. 

For illustrative purposes, we consider a spherically symmetric baryonless 
drop of quark-gluon plasma. The main conclusions derived from our anal- 
ysis would neither depend on the exact shape of a QGP drop nor on the 
details of the equation of state such as nonzero baryon density. Instead, it 
is the presence of the boundary, which is the origin and the reason of the 
distinctive properties of our model. 

So we consider a spherical drop at a constant temperature Tq and pressure 
po in the initial moment. The initial velocity is assumed to rise linearly 
from |v| = at the center of the drop to |v| = vq at the boundary. We have 
tried widely different values of po and vq. We have also tested the sensitivity 
of the results on the choice of the bag constant B in the equations of state 
and to the strength of the surface tension, which determines the mechanical 
stability of the drop against breaking. 

Consider first the case of zero initial velocity vq = (fig. 1) and following 
parameters: radius rg = 3 fm, po = 6.5 GeV/ fm^, B = 0.2 GeV/ frrfi. 
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The expansion of a QGP drop starts from the nearby-surface regions. The 
overfall of pressure (from plasma to vacuum) leads to the formation of 
a rarefaction wave, which propagates from the surface of the drop to its 
center (fig. 1 a,b,c). The velocity of this wave is the speed of sound. This 
wave puts the plasma in motion, thus causing the process of expansion. 
The continuous expansion of the QGP drop during all this time results in 
a monotonic decrease of the pressure and the temperature in the whole of 
the drop (fig. 1 a-h). 

If the initial velocity is high (vq > Vg, Vg being the speed of sound) the 
evolution starts from the creation of excess pressure region at the boundary 
of the QGP drop. This leads to a formation of a compression wave, which 
also moves from the surface of the drop to its center with the speed of 
sound. In the same time the QGP drop expands as a whole, so that the 
pressure and the temperature decrease as in the previous case. When 
the compression wave reaches the center of the drop and reflects from it, 
the evolution of the drop becomes very similar to the case of zero initial 
velocity. 

The expansion of a QGP drop may be induced either by the rarefaction 
wave (if vq = 0) or the initial motion (vq > 0), or both these reasons, but 
the behaviour of plasma looks the same in the sense that the expansion 
and cooling are always accompanied by the propagation of waves. 

Depending on the initial amount of energy (or the initial temperature) 
the expansion may take different time. The compression waves may pass 
through the drop several times, reflecting successively from the center and 
from the boundary. The expansion and cooling of plasma carry on until 
the pressure falls below zero or the conditions for a phase transition are 
reached. What happens next depends on the mechanical and thermody- 
namical stability of the plasma. For the time, let us forget about phase 
transitions and concentrate on a purely hydrodynamical problem. 

The mechanical stability of plasma depends on the strength of its surface 
tension. The estimations of this parameter are uncertain, and hence we 
will consider two opposite extreme cases. 

First, assume that the surface tension is always strong enough to prevent 
the QGP drop from breaking into smaller droplets. Then the expansion of 
plasma by inertia may lead to states with negative pressure. If the initial 
pressure was small compared to the bag constant {po < B/2), the negative 
pressure appears for the first time at the moment when the rarefaction 
wave reaches the center of the QGP drop. Since that moment the region of 
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negative pressure spreads with the reflected wave towards the drop bound- 
ary. The moment when it reaches the surface of the drop corresponds to 
the maximal expansion of plasma and to its deepest overcooling. At this 
time the pressure is almost uniformly negative in all the volume of the drop 
and its value is approximately equal to the initial pressure taken with the 
opposite sign: pq ~ —pQ. After that the drop starts to squeeze. At the 
end of the period of squeezing plasma returns to the state, which is similar 
to the initial one, both in its pressure (or temperature) and the occupied 
volume. 

If the initial pressure is high (po > B/2), the process of expansion and 
cooling takes a longer time. As result of reflaction of a plasma flow from 
the drop boundary, some excess of pressure is formed near the boundary. 
It creates the compression wave passing through the plasma drop many 
times reflecting successively from its center and from its surface at the 
background of the more strong rarefaction process. This compression wave 
slows down the expansion of plasma and eventually stops it. It is remark- 
able that even at the moment of maximal overcooling, the negative pressure 
can never fall below a definite value, which depends on the bag constant B 
but does not depend on the initial pressure po: p > Pmin, Pmin — —B/2 
(fig. 1 j). The subsequent squeezing of plasma restores positive pressure 
in the whole of the drop. However, the value of this positive pressure pg 
is not equal to the initial value j^o, but is nearly the minimal (negative) 
pressure taken with the opposite sign: Pq ~ —Pmin (fig- 1 !)• Neither does 
the occupied volume contract to the initial size. 

If to allow a further evolution of plasma, one would see a second period 
of expanding and squeezing. It essentially resembles the first period, but 
corresponds to a reduced initial pressure, i.e. p'q instead of po (figs. 11). At 
the end of this period plasma comes to a state with Pq — p'q, as is typical for 
processes starting from relatively low pressure. All the subsequent periods 
are simply a repetition of this picture. 

On the other hand, if the plasma surface tension is weak the evolution 
of the drop comes to end as soon as the pressure falls to zero. Then the 
creation of vacuum bubbles (or breaks) comes into play. The breaks prevent 
the temperature and the pressure from further decreasing and "freeze" the 
state with ]7 ~ 0. 

Most probably, the vacuum bubbles will break the QGP drop into smaller 
droplets, but their summary volume must conserve. The volume cannot 
increase because this would lead to further decreasing pressure and will 
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simply end up with new breaks. The volume cannot decrease because 
there is no force, which could cause squeezing of plasma, while the motion 
by inertia tends to move the drops away from each other. Hence, the 
primordial QGP drop will transform into a sponge- like object or into a 
cloud of smaller droplets. The pressure in this system freezes at the ^ ~ 
level. 

It is important to point out that the dynamical limit for overcooling 
Pmin — —B/2 remains valid in all the cases. The existence of this limit 
restricts the possible ways of hadronization. An explanatory overview of 
all types of phase transitions may be found in the textbook [jT^], and so 
we will not repeat it here. According to the estimations made in jS], 0] the 
supersonic condensation and the detonation may only occur if p < —B /2. 
Since this degree of the overcooling can never be achieved, no matter what 
kind of initial conditions are chosen, we are left with slow processes like 
equilibrium phase transition and deflagration. The deflagration must be 
considered the most realistic scenario for the hadronization of plasma. 

Now let us take the deflagration into account. By the definition [jT^], the 
deflagration is a slow process, which spreads with a velocity much less than 
the speed of sound. There are no reasons to accept that the plasma can 
convert into hadrons in a moment, once the pressure has fallen down to 
the equilibrium threshold. Moreover, the deflagration takes place only at 
the surface of plasma drop, but not in its inner regions. 

In the present consideration we have simulated the deflagration as a 
process, which starts at the outer surface of the drop and spreads towards 
the center with a constant velocity Vdef = 0.1 c Q. The deflagration 
proceeds during the time when the necessary condition < j9 < pth is 
satisfied, where the threshold value pth has been estimated in [§,0]: pth = 
30 MeV/fml 

Fig. 2 illustrates the evolution of a QGP drop with the account of de- 
flagration and ^0 = 2.56GeV/ fm^^ (assuming the negative pressure is 
allowed). Hadronization reduces the mass of plasma, but proceeds with a 
relatively small velocity. One can see several oscillations of the drop be- 
fore the termination of the phase transformation. If the initial drop would 
break into smaller droplets (assuming the negative pressure is unstable) 
the deflagration would proceed faster because of enlarged total surface of 
the system. However, it takes a rather long time to reach even the state, 
where the deflagration is possible. This fact has a direct consequence in 
the enhanced yield of thermal photons and lepton pairs. 
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5. The yield of thermal photons 

To estimate the yield of photons from QGP we recall the analytical re- 
sults of ref. [TTj. The emission of photons is mainly due to two elemen- 
tary processes, i.e. the gluon Compton scattering and the annihilation of 
quarks: 

q{Pi)+9{P2)^l{k) + q{k'), (12) 
q{pi) + q{V2)^l{k)+g{k'). (13) 

Here the letters in the parentheses denote the momenta of the particles. 
In the lowest QCD order, the spin- and color-averaged matrix elements for 
these processes read: 



IGtt^ 
1287r2 




W = 7^--^ + 77— I (15) 

where Qf and m are the electric charge and the mass of a quark (/ = 
s), and s, t and u are the standard Mandelstamm variables: s = 

{pi + ^32)^ t = (pi- kf, u = {pi- Jzf. 

The total number of photons produced in unit volume of plasma per unit 
time is given by: 

'^1^ = (2^/2E^"^^^/2Er"^^^V^-(^^)x 

S\pi+p2-k-k') j:\Tabf (16) 

a, 6 

where cu is the photon energy, £"1, £"2, E' are the energies of the particles 
with momenta pi, ^2, p', respectively; A^^ is the number of degrees of free- 
dom (spin, color) in the initial state of the reaction, n{E) stand for Fermi 
or Bose distributions, and n{E) = 1 — n{E). If to use an approximate 
Maxwell-Boltzmann formula n{E) = exp(— £"/£), h{E) = 1 instead of ex- 
act Fermi and Bose distributions, the multidimensional integration in ([T^) 
may be performed analytically. The yield of photons per unit volume and 
unit time is then 

TT 9T 

N = -{T.Q})c^c.s{T)H—)T\ (17) 
6 f m 
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Here the running QCD coupling constant is 

^s{T) = — ^f^. rp/. (18) 

(33 - nf)ln{KT/AMs) 

According to lattice calculations, ~ 4 and Tc/ Ams — 1-8. 

Equation (jlBD leads to the following pseudorapidity and the transverse 
momentum distributions of photons: 

H"^ N 4 Ak T h 

^ = ^(E Q?) « am r- ln(-^) exp(-^) cosh(, - ,')■ (19) 

To calculate the total yield of photons from plasma we first integrate the 
hydrodynamical equations of motion and find the temperature at every 
time at every space point. Then we use eqs. (|r7]) and 

For comparison, the prediction of Bjorken model for the evolution of 
temperature with time is 

T = To{To/Tr% (20) 

where n = 1, 2, 3 refer to one-, two- or three-dimensional expansion of 
plasma. 

The results of our calculations are seen in fig. 3. Fig. 3a shows the 
total rate of photon emission as a function of time. Figs. 3b, c show the 
distribution of the produced photons on the transverse momentum for = 
258MeV and To = 319MeV respectively. Solid curves represent integrated 
spectra from the beginning till the moment when the pressure reaches zero 
level in the whole volume of the drop. Dashed curves correspond to the 
emission of photons in accord to Bjorken scenario. Dotted curves refer to 
photons originating from decays of final-state hadrons, as is estimated in 

Variations in the model parameters may shift the presented numbers but 
cannot smash the general picture. The conclusion on the long lifetime of 
plasma and on the enhanced yield of photons remains valid within a wide 
range of model assumptions. 



6. Conclusions 

We have established a novel mechanism of expansion and cooling of 
QGP, which consists in multiple propagation of rarefaction and compres- 
sion waves trough the plasma. Since the speed of these waves is almost 
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constant and equal to the speed of sound, the time needed to cool plasma 
increases with the size of QGP drop. 

Even if to allow the existence of states with negative pressure, the hydro- 
dynamical evolution of plasma cannot lead the pressure below a definite 
limit, which is determined by the bag constant in the equations of state: 

Pmin — B /2. 

Since this degree of overcooling is insufficient for supersonic phase transi- 
tions, hadronization may only occur due to rather low processes like defla- 
gration or equilibrium transition. As a consequence, the lifetime of plasma 
extends considerably. 

Both the slow cooling mechanism and slow hadronization process result 
in a longer period of photon radiation. The total yield of thermal photons 
gets comparable with or even exceeds the contribution from hadron decays. 

The enhanced yield of photons is a fundamental result, which may be 
proven by a direct experimental measurement. The impossibility of su- 
personic phase transitions from plasma to hadrons puts some restrictions 
on the proportion between different hadrons and, hopefully, may also be 
experimentally tested. 
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Figure captions 
Fig. 1 

Pressure as a function of the coordinate r along the radius of a QGP drop, 
plotted for different times (t). The initial radius of the drop, the pres- 
sure and the velocity are: ro = 3 fm, po = 6.5 GeV/fm^, vq = and 
B = 0.2 GeV/fm^. Mind the different ordinate scale in different plots 
(a-1). 

Fig. 2 

The process of expansion and cooling of a QGP drop with the account of 
deflagration. 

Fig. 3 

a - Total rate of photon emission, Tq = 258 MeV. 
b - Transverse momentum distribution of photons. 
Solid curve - calculations within the present model. 
Dashed curve - calculations according to Bjorken model. 
Dotted curve - background from decays of hadrons. 
The initial temperature Tq = 258 MeV. 
c - The same as in b, but for Tq = 319 MeV. 
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